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Abstract 

Motivated by recent Pierre Auger result on the correlation of the highest-energy cosmic rays with 
the nearby active galactic nuclei, we explore possible ultrahigh energy cosmic ray (UHECR) source 
distributions and their effects on GZK horizons. Effects on GZK horizons by local over-density of 
UHECR sources are examined carefully with constraints on the degree of local over-density inferred 
from the measured UHECR spectrum. We include the energy calibration effect on the Pierre Auger 
data in our studies. We propose possible local over-densities of UHECR sources which are testable 
in the future cosmic ray astronomy. 

PACS numbers: 95.85.Ry, 96.50.sb, 96.50.Vg 
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I. INTRODUCTION 



Recently, Pierre Auger observatory published results on correlation of the highest-energy 
cosmic rays with the positions of nearby active galactic nuclei (AGN) . Such a corre- 
ation is confirmed by the data of Yakutsk while it is not found in the analysis by HiRes 
J]. In the Auger result, the correlation is maximal for the threshold energy of cosmic rays 
at 5.7 X 10^^ eV, the maximal distance of AGN at 71 Mpc and the maximal angular sepa- 
ration of cosmic ray events at ip = 3.2°. With the same threshold energy, and the angular 
separation < Q°, the correlation remains strong for a range of maximal AGN distance 
between 50 Mpc and 100 Mpc. Due to increasing efforts on verifying the Auger result, it is 
worthwhile to examine the above correlation from a phenomenological point of view. 

Since the angular scale of the observed correlation is a few degrees, one expects that these 
cosmic ray particles arepredominantly light nuclei. The effect of GZK attenuations on these 
cosmic ray particles [s], [g] can be described by a distance scale referred to as "GZK horizon" 
which is a function of the selected energy threshold for the arriving cosmic ray particles. By 
definition, the GZK horizon associated with a threshold energy Et^ is the radius of a spherical 
region which is centered at the Earth and produce 90% of UHECR events arriving on Earth 
with energies above Eth_- With continuous energy loss approximation, the GZK horizon for 
protons with E'th = 57 EeV is about 200 Mpc by assuming a uniformly distributed UHECR 
sources with identical cosmic ray luminosity and spectral index The calculations based 
upon kinetic equation approach or stochastic energy loss also reach to similar conclusions 

,3- 
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The departure o 
of the V-C catalog 



theoretically calculated GZK horizon to the maximum valid distance 



lOj employed in Pierre Auger's analysis, which is around 100 Mpc, can 



be attributed to several factors. As mentioned in |2i], such a deviation may arise from non- 
uniformities of spatial distribution, intrinsic luminosity and spectral index of local AGN. In 
addition, the energy calibration also plays a crucial role since the GZK horizon is highly 
sensitive to the threshold energy Eth- Energy values corresponding to the dip and the GZK 



cutoff of U 
periments 



3ECR spectrum were used to calibrate energy scales of different cosmic ray ex- 
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12| . It has been shown that all measured UHECR energy spectra can be 



wrought into good agreements by suitably adjusting the energy scale of each experiment 
ll| . Keeping the HiRes energy scale unchanged, the energy-adjustment factor A is found to 
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be 1.2, 0.75, 0.83 and 0.625 respectively for Auger, AGASA, Akeno and Yakutsk. Further- 
more, it has been shown that a different shower energy reconstruction method infers a 30% 
higher UHECR energy than that determined by Auger's fluorescence detector-based shower 
reconstruction 1^. 

In this paper, we investigate the consistency between Auger's UHECR correlation study 
and its spectrum measurement. As just stated, the V-C catalog used by Pierre Auger for 
the correlation study is complete only up to 100 Mpc while the GZK horizon for Eth = 57 
EeV is generally of the order 200 Mpc. We first consider the local over-density of UHECR 
sources as a possible resolution to the above discrepancy. It is motivated by the existence of 
Local Supercluster (LS) which has a diameter of the order 60 Mpc. In LS, the over-density 



1^. 
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17( 1 to account for 



of galaxies has been estimated to be ~ 2 

The local over-density of UHECR sources has been invoked 
AGASA data 18|, ll9|. Such a density distribution naturally leads to a smaller GZK horizon. 
However, it also significantly affects the UHECR energy s pec trum in (5 ■ 10^^ — 10^°) eV 
region. Hence fittings to the measured UHECR spectrum 20|] can provide information on 
the degree of local over-density. Subsequently, the magnitude of GZK horizon can be better 
estimated. 

We next study the energy calibration effect on the estimation of GZK horizon and the 
spectrum of UHECR. Certainly a 20% — 30% upward shift on UHECR energies reduces the 
departure of theoretically calculated GZK horizon to the maximum valid distance of V-C 
catalog [2I. The further implications of this shift will be studied in fittings to the shifted 
Auger spectrum. 

We fit the UHECR spectrum for events with energies above 10^^ eV. This is the energy 
range where the GZK attenuation exhibits its effect. It is also the energy range where the 
local over-density of UHECR sources shows significant effects. In our analysis, we take the 
UHECR as protons, which is hinted in the Auger events with energies > 57 EeV although 



;he composition study by the same group suggests a heavier composition 



br E < 40 EeV 



2l| . The HiRes experiment measures the composition up to 50 EeV [22i] and obtains a 
composition lighter than that of Auger. For E > 50 EeV, the event number is still too 
small for the composition study. To fit the UHECR spectrum at the highest energy, it is 
more appropriate to treat the cosmic ray energy loss as a stochastic process [23|. There 
are numerical packages available for treating stochastic energy loss of cosmic ray particles 
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24 |25|. We employ the latter package for our calculations. Although UHECR loses its 



energy mostly by scattering off CMB photons, it also 
scattering off infrared background photons 
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OSes some amount of energy by 



301 1 . Thus we include the infrared 
photon contribution to the UHECR energy attenuation. Source evolution n{z) = nQ{l + z)^ 
is adopted in the calculation of GZK horizon and spectrum, where hq is the source number 
density at the present epoch. It is from the generally-accepted soft evolution model which 



3l|. 



traces the star formation history and has been adopted in previous works 

We discuss about GZK horizons in Sec. II. We calculate the accumulated event probabil- 
ities of UHECR for Eth = 57 EeV, 70 EeV, 80 EeV and 90 EeV respectively. GZK horizons 
corresponding to different i?th are tabulated. We also calculate GZK horizons with local over- 
density of UHECR sources taken into account. In Sec. HI, we fit the measured UHECR 
spectrum with various local over-densities of UHECR sources and obtain information on the 
degree of local over-density. To study the energy calibration effect, we also perform fittings 
to the shifted UHECR spectrum. Sec. IV contains discussions and conclusions. 

II. THE ACCUMULATIVE EVENT PROBABILITIES OF UHECR 



(1) 



For single UHECR source, the cosmic-ray energy attenuation is governed by the equation 

dME,t) _ d 
dt dE 

in the continuous energy loss approximation. This equation results from the number con- 
servation of cosmic-ray particles in the energy attenuation process. The cosmic-ray energy 
loss per unit time —dE/dt is due to the cosmic expansion and its scattering with cosmic 
microwave background photons through photo-pion production process P7 — > Nn and pair 
production process P7 Pe^e^. The above attenuation equation is well known [32]. In the 
current context, the solution of Eq. ([1]) can be expressed in terms of the red-shift variable 
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{E, Zs) X exp 



dz' 



[1 + z') dboiil + z')E) , 1 

X , ^ — - — - + 



(2) 



H{z') dE ' l + z\ 

where Zg is the red-shift of the UHECR source and the function 60 is related to the rate of 
cosmic-ray energy loss at the present epoch by 

dE 



dt 



{z = 0) = bo{E) + HoE, 



(3) 



where Hq is the present value of Hubble constant. The UHECR has an energy E at the 
source with red-shift Zs and its energy is downgraded to E at the red-shift z. The energy E 
is a function of E and z so that E{E, Zs) = E and 

E 



dE _ bo ((1 + z)E) 



dz 



Hz) 



l + z 



(4) 



Due to the non-trivial form of bo, one resorts to numerical methods for computing the 
function E and the flux (j)i^{E,z). 

We have mentioned that the stochastic nature of UHECR energy loss can not be over- 
looked for shorter propagation distances One then treats the energy attenuation by 
photo-pion production as a stochastic process while treating other attenuations as continu- 
ous processes. 



Eth=57EeV 



Ech=7 0EeV 





Eth^BOEeV 



Eti,= 90EeV 





FIG. 1: The accumulative event probability P{D, E'th) as a function of D for E'th = 57 EeV, 70 EeV, 
80 EeV and 90 EeV respectively. The horizontal dash line in each panel denotes P{D, -Eth) = 0.9. 
The red, green, blue and black curves represent results from models with over-density n{l < 
30Mpc)/no = 1, 2, 4, and 10 respectively. The intrinsic spectrum index 7 = 2.4, energy cut 
-Scut = 1000 EeV and the source evolution model n{z) = no(l + -z)^ are used for calculations. 

To facilitate our discussions, we define the accumulative event probabilities of UHECR 
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as 



dl ■ N{1, E, 



P(D,E, 



th 



thj 



j^^di-N{i,E,^y 

where N{1, Et-^) ■ dl is the number of cosmic ray events which are originated from sources 
at distances between I and / + dl from the Earth and arrive at the detector with energies 
above Eth- We calculate P{D, Eth) for various local over-densities of UHECR sources. The 
source distribution over the red-shift is taken as n{z) =^0(1 + 2;)^ and the energy spectrum 
of each source is taken to be the form, (p^iE) = dN/dE = AE'"', with the maximal energy 
Ecnt = 1000 EeV. We choose 7 = 2.4, 2.5 and 2.6 where 7 = 2.5 gives the best fitting to 
the measured UHECR spectrum as will be shown in the next section. The accumulative 
event probability P{D, Eth) for Eth = 57 EeV, 70 EeV, 80 EeV and 90 EeV are shown in 
Fig. [T]for 7 = 2.4. Results for 7 = 2.5 and 7 = 2.6 are not distinguishable from those for 
7 = 2.4. In each panel, the red, green, blue, and black curves represent local over-density 
n(/ < 30Mpc)/no = 1, 2, 4, and 10 respectively. The local over-density n{l < 30Mpc)/no = 
k is defined explicitly as 



n(/ < 30Mpc)/no = k{l + zf, 
n{l > 30Mpc)/no = (1 + zf. 



(6) 



The horizontal dash line in each panel denotes P{D, Eth) = 0.9. The intersection of this line 
with each color curve gives the GZK horizon corresponding to a specific local over-density 
characterized by the ratio n{l < 30Mpc)/no. 



TABLE I: GZK horizons of UHECR calculated with the local over-density n{l < 30Mpc)/no = 
1, 2, 4, and 10, and arrival threshold energy Eth = 57 EeV, 70 EeV, 80 EeV and 90 EeV respec- 
tively. The listed numbers are in units of Mpc. 



n{l < 30Mpc)/no 


Eth = 57 EeV 


^th = 70 EeV 


Eth = 80 EeV 


Eth = 90 EeV 


1 


220 


150 


115 


90 


2 


210 


140 


105 


75 


4 


195 


120 


85 


60 


10 


155 


85 


50 


30 



GZK horizons corresponding to different local over-densities and Eth are summarized in 
Table I. It is seen that local over-densities up to n{l < 30Mpc)/no = 4 do not alter GZK 
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horizons significantly for a given Eth- One could consider possibilities for higher local over- 
densities. However, there are no evidences for such over- densities either from astronomical 
observations Q or fro,„ fltturgs to the .measured UHECR speetrtun. We .rote that GZK 

horizons are rather sensitive to ii^th- Table I shows that GZK horizons are ~ 100 Mpc or 
less for > 80 EeV. 



III. FITTINGS TO THE UHECR SPECTRUM MEASURED BY PIERRE AUGER 



As mentioned earlier, the local over-density of UHECR sources affects the cosmic-ray 
spectrum at the highest energy, especially at energies higher than 5 ■ 10^^ eV. Hence the 
degree of local over-density can be examined through fittings to the measured UHECR 
spectrum as will be shown momentarily. 

Fittings to the Auger spectrum have been performed in 
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34j . l35l . l36l|.. In our work. 



we take into account the over-density of UHECR sources in the distance scale / < 30 
Mpc. As stated previously, we take the UHECR to be protons. Figure [2] shows our 



TABLE II: The values of total from fittings to the Auger measured UHECR spectrum. Numbers 
in the parenthesis are values from fittings to the 8 data points in the energy range 19.05 < 
logiQ{E/eV) < 19.75. The last 4 data points record events with energy greater than 71 EeV. 



n{l < 30Mpc)/no 


1 


2 


4 


10 


7 = 2.5 


14.12(9.34) 


14.61(9.93) 


17.09(10.50) 


28.09(13.93) 


7 = 2.6 


16.64(12.28) 


15.56(11.90) 


16.01(11.83) 


20.76(11.67) 



fittings to the Auger measured UHECR spectrum with 7 = 2.5 and 2.6 respectively. We 
take the red-shift dependence of the source density as n{z) = no{l + z)"^ with m = 3. 
We have fitted 12 Auger data points beginning at the energy 10^^ eV. We make a fiux 
normalization at 10^^ eV while varying the power index 7 and the the degree of local over- 
density, n{l < 30Mpc)/no. Part of values from our fittings are summarized in Table II. We 
found that 7 = 2.5, n{l < 30Mpc)/no = 1 gives the smallest value with x^/d.o.f. = 1.57. 
For the same power index, the large local over-density n{l < 30Mpc)/no = 10 is ruled out 
at the significance level a = 0.001. For 7 = 2.6, n{l < 30Mpc)/no = 10 is ruled out at the 
significance level a = 0.02. 
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FIG. 2: Fittings to the Auger measured UHECR spectrum where the red, green, blue and black 
curves denote the model with the local over-density n(/ < 30Mpc)/no = 1, 2, 4, and 10 respectively. 
Solid curves correspond to 7 = 2.6 while dash curves correspond to 7 = 2.5. We take the source 
evolution parameter m = 3 throughout the calculations. 

We note that, for both 7 = 2.5 and 7 = 2.6, the GZK horizon with n{l < 30Mpc) /% = 10, 
Etk = 57 EeV, m = 3 and E^^t = 1000 EeV is about 155 Mpc. Since n{l < 30Mpc)/no = 10 
is clearly disfavored by the spectrum fitting, one expects a GZK horizon significantly larger 
than 155 Mpc for Eth = 57 EeV. 

We next perform fittings to the shifted Auger spectrum. The results are shown in Fig. [3] 
where the cosmic ray energy is shifted upward by 30%. Part of values are summarized in 
Table III. The smallest value occurs approximately at 7 = 2.4, n{l < 30Mpc)/no = 2 with 
xVd-o.f = 0.82. For 7 = 2.5, xVd-o.f = 1.31, 0.96 and 0.87 for n{l < 30Mpc)/no = 1, 2 and 
4 respectively. It is seen that values from current fittings are considerably smaller than 
those from fittings to the unshifted spectrum. Given a significance level a = 0.1, it is seen 
that every local over-density listed in Table III except n{l < 30Mpc)/r;,o = 10 is consistent 
with the measured UHECR spectrum. It is intriguing to test such local over-densities as 
will be discussed in the next section. 

We note that, with a 30% upward shift of energies, the cosmic ray events analyzed in 
Auger's correlation study would have energies higher than 74 EeV instead of 57 EeV. The 
GZK horizon corresponding to Et^^ = 74 EeV is 120 Mpc for n{l < 30Mpc)/no = 2 and 105 
Mpc for n{l < 30Mpc) /no = 4. 

We have so far confined our discussions at m = 3. In the literature, m has been taken as 
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Log 10 [Energy (eV) ] 

FIG. 3: Fittings to the Auger measured UHECR spectrum with a 30% upward shift on UHECR 
energies where the red, green, blue and black curves denote the model with the local over-density 
n{l < 30Mpc)/no = 1, 2, 4, and 10 respectively. Solid curves correspond to 7 = 2.4 while dash 
curves correspond to 7 = 2.5. We take the source evolution parameter m = 3 throughout the 
calculations. 

TABLE III: The total values from fittings to the Auger measured UHECR spectrum with a 
30% upward shift on UHECR energies. Numbers in the parenthesis are values from fittings to 
the 8 data points in the energy range 19.16 < log]^Q(£'/eV) < 19.86. The last 4 data points record 
events with energy greater than 92 EeV. 



n{l < 30Mpc)/no 


1 


2 


4 


10 


7 = 2.4 


8.65(4.30) 


7.39(4.67) 


10.26(6.35) 


27.31(13.34) 


7 = 2.5 


11.82(6.16) 


8.67(5.49) 


7.78(5.23) 


16.18(7.39) 



any number between and 5. It is demonstrated that the effect on UHECR spectrum caused 
by varying m can be compensated by suitably adjusting the power index 7 3l||- Since GZK 
horizons are not sensitive to 7 and m, results from the above analysis also hold for other 



m s. 
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IV. DISCUSSIONS AND CONCLUSIONS 



We have investigated the consistency between Auger's latest result on the correlation 
of UHECR sources with positions of nearby extra-galactic AGN and its measured UHECR 
spectrum. As stated before, this investigation is motivated by the fact that the V-C catalog 
used by Pierre Auger for the correlation study is reliable only up to 100 Mpc while the 
GZK horizon for Et^ = 57 EeV is generally of the order 200 Mpc. We have explored the 
possibility for local over-density of UHECR sources, which is expected to shorten the GZK 
horizon for a given threshold energy of arrival cosmic-ray particles. This is indeed the case 
as can be seen from Table I. On the other hand, the effect is far from sufficient to shorten 
the GZK horizon at E'th = 57 EeV to ~ 100 Mpc for a local over-density of UHECR sources 
consistent with the measured UHECR spectrum. 

We have performed a upward energy shift to the Auger measured UHECR spectrum. As 
said, a upward energy shift is motivated by simulations of shower energy reconstructions as 
well as the requirement of reproducing the theoretically predicted GZK cutoff energy. With 
a 30% energy shift, each cosmic ray event used by Auger for the correlation study would 
have an energy above 74 EeV instead of 57 EeV. GZK horizons corresponding to E'th = 74 
EeV then match well with the maximum valid distance of V-C catalog. Fittings to the 
shifted Auger spectrum indicate a possibility for the local over-density of UHECR sources. 

We point out that the local over-density of UHECR sources is testable in the future 
cosmic ray astronomy where directions and distances of UHECR sources can be determined. 
Table IV shows percentages of cosmic ray events that come from sources within 30 Mpc for 
different values of E'th and n{l < 30Mpc)/no. We take 7 = 2.4, m — 3 and E^ut — 1000 
EeV for calculating these percentages. We note that these percentages are not sensitive to 
the above parameters. For i?th = 57 EeV, only 17% of cosmic ray events come from sources 
less than 30 Mpc away for n(/ < 30Mpc)/no = 1. For n{l < 30Mpc)/no = 2 and the same 
threshold energy, 30% of cosmic ray events are originated from sources in the same region. 

It should be stressed that we have focused only on resolving the apparent discrepancy 
between the GZK horizon at E^^ — 57 EeV and the maximum vahd distance of V-C catalog. 
The statistics analysis for establishing the source correlation is an independent issue beyond 
the scope of the current paper. We have found that the above discrepancy can not be 
resolved by merely introducing the local over-density of UHECR sources. On the other hand. 
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TABLE IV: Percentages of cosmic ray events that come from sources within 30 Mpc for different 
values of Eti^ and local over-density n{l < 30Mpc)/no. 



n{l < 30Mpc)/no 


Eti, = 57EeV 


= 70 EeV 


Etk = 80 EeV 


Eti, = 90 EeV 


1 


0.17 


0.27 


0.36 


0.46 


2 


0.30 


0.43 


0.53 


0.63 


4 


0.46 


0.60 


0.70 


0.77 


10 


0.68 


0.79 


0.85 


0.89 



if Auger's energy calibration indeed underestimates the UHECR energy, such a discrepancy 
can be reduced. More importantly, fittings to the shifted Auger spectrum indicate a possible 
local over-density of UHECR sources, which is testable in the future cosmic ray astronomy. 
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